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On the Theory of Stark Broadening of Spectral Lines in a Plasma* 

G. K. Oertel 

MSA, Langley Research Center 
ABSTRACT 

The theory of Stark broadening of isolated spectral lines by 
electrons in a plasma was developed by Griem et al (ref. l), and numer- 
ical results were obtained for a large number of neutral and singly 
ionized spectral lines (ref. 2). Large discrepancies between measured 
and predicted widths of ion lines (e.g., Jalufka et al, ref. 5 and 4) 
have led to a re investigation of the underlying assumptions. It is shown 
that the calculated widths do not change by more than a few percent as a 
result of certain corrections. The effect of different strong collision 
cutoff procedures is as important but still far too small to explain the 
experimental results for ions. Debye shielding must often be taken into 
account in the electron density determinations from broadening of neutral 
lines in typical experiments and was therefore built into a conputer 
program for electron impact widths. The resiilts are compared with estimates 
on the basis of the usual procedure (ref. l) for estimating Debye shielding 
effects. 

An abstract of this paper was presented at 8th Annuial Meeting of the 
American Physical Soc., Div. of Plasma Pl:orsics; Boston, Ma.ss.) Nov. 2-^, 1966. 
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I. lOTRODUOPIOK 

A sizable discrepancy between experiment and theory in the Stark 
broadening of lines emitted by sing3y charged argon ions has been 
reported previously (refs. 5 and 4 ). The measured lines were wider by 
factors of typically 2.5« 

The Stark broadening theory for isolated lines has been reviewed 
critically in an effort to determine the cause of this discrepancy. 

J. Cooper^ and the author are cooperating in the development of a gen- 
eralization of the relevant Stark broadening theory (ref . 1) which has 
succeeded in removing the discrepancy for all cases for which numerical 
results are now available (for a different approach see also ref. 5)» 
While it would be premature to present details of this most recent 
development it is considered to be appropriate to report some refinements 
to the theory as published in reference 1 (GBKD). 


II. THE STARK BROADENING FUNCTIONS 

The functions A(z) and B(z) which arise from the first 
nonvanishing term in the Djrson series in GBKD are defined by 


A(z) + iB(z) = 




dx. 


dxg e 


iz(xi-X2) 


f(x^,X2) (1) 


^at JILA, Boulder, Colorado 
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irith 


f(X]^,X2> « I 

Width and shift due to electron iu^pact broadening are proportional to 
thermal averages over the functions a(z) and b(z), respectively, which 
are defined by 



WinPH 



( 2 ) 



(3) 


While a(z) had been expressed in terms of modified Bessel functions 
in GBKO the in^iortant function a(z) had been obtained numerically on a 
digital computer. It txjrns out that the integral can be evaluated 
analytically, and that the resulting a(z) is consistently smaller than 
used in GBKDj 

a(z) = z Kq(|z|) Kj^(Iz I) (4) 

As in GBKO, K^ and are modified Bessel functions. Figure 1 shows 
a(z) along with its asymptotic forms for small and large z and with the 
numerical results of GBKO for comparison. 

Since most of the broadening comes from z « 1 where the correct 
asymptotic form was used in GBKD one would expect little effect on the 
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nTmibers coiapxibed there and in later reports by Griem. This is illustrated 
hy table 1 which shows some widths computed with the old and new functions. 
The change of abottt five percent or less is probably negligible compared 
to the uncertainties in the matrix elements. It may also be noted that a 
comparable uncertainty is introduced by the choice of the strong collision 
cutoff. 

TABLE 1 



mums [a] AT 10 ^^ ELECTROHS/CM^ 


a(z) 

Promj 

GBKO^ 

Present Form 

(cutoff I (see footnote 2) 

.17, 

,l 6 u 

HE 5876 < 

II (see footnote 3) 



1 " 

.18- 

3 


i " 

I 

•360 

•5^8 

HE 5016 < 


1 " 

II 





5 

i " 

I 

.OID. 

.009- 

AR II < 


0 

5 

V" 

II 

.0122 

.Ollg 


The function B(z) is obtained numerically via the principal value 
integral 

1 

as recalculated by this author. 

2 

as in reference 1. 

5as in reference 2. 
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B(z) 


2z 





(5) 


Revised values of B(z) between 0 ^ z ^ 5 are within typically one 
percent of those in GBKD. However, b(z) had to be Increased by nearly 
20 percent for small values of z. This Is illustrated in figure 2 where 
B(z) and b(z) are plotted along with b(z) as given in GBKD. 

Unlike the values in GBKD the revised values were obtained by inte- 
grating over from z to infinity. An excellent check is provided 

by the known value b(o) s | (GBKD). The revised value is equal to jt/2 
within better than one percent although (see eq. 5) the integration over 
B(z)/z must be carried out over the entire range of z (0 to «) to obtain 
b(o). Although the new values differ most from the GBBD values for z < 2, 
the resulting changes by typically 10 percent in the shift are not critical 
because of the inherently much larger uncertainties in all confuted shifts. 

The successful redetermination of b(z) was made possible by a much 
extended asymptotic e3q>ansion of B(z) and b(z). It can be shown that 
the asymptotic e:^)ansions of B(z) and b(z) for large z are 


B(z) 


k 

E 

n=l 


a 

n 

,2n-l 


( 6 ) 


and 


b(z) = 


k 

E 


a. 


n 


(2n - l)z 


2n-l 


(7) 
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respectively, with 


-i: 


x^,Xg 


dx. 


2n-2 


( 8 ) 


^2"^1 


f ^x^,Xg^ is the same as in equation 1. 

It can now he shown with a considerable amount of algebra that 


jt 


n-1 nH 

2 2 

n (2i+l) Z 
i=l 


2 \ i-1 n-i+2 / T V 


vi-li 


n+1 


- i 


i-1 


X n (nt2j) n (nH-5-2j)| 
^=1 J =1 


ntl n»l /n-l^ 

2 . 2/2 

n (2i-l)(-i)'- Z 

i=l i=l \i-l. 


n+1 


- 1 


i-1 


X n (n+2j) n (n+l-2j) 
J=1 j=l 


n+l-i 


n+2-i 

" 

N 

" ( 

1 - 

1- n C 

1 - i-1 


j=i ^ 

2j; 


2jy/ 



( 9 ) 


This expression is best evaluated on a conputer. The first few coefficients 
are listed in table 2. 
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0?ABLE 2 


CX)EFPICIEI®?S OP THE ASXMPTCTIC EXPAESION OP B(z) 


n 


n 


1 

3 

5 

7 

9 


jt/^ = •T85^ 

9it/52 » .8856 

42.28 

750.4 


At z = 5 the asymptotic value for B(z) to order z ^ agrees with 
the revised value frcm the principle value integration within .02 percent. 

III. DEBYE SHIELDING 


Debye shielding is often only a correction and may be taken into 

account approximately by the procedure outlined by GBKD. A shielded width 

w may be estimated from the tabulated unshielded width w by 
s o 



where co , is the separation (in angular frequency tinlts) of the closest 
cax 

perturbing level, the plasma frequency, v the mean thermal electron 
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speed, and the minimum iuqpaet parameter. Sim;e j should he 

( max \ 

z I J , vhich means that the integration over impact parameters 
p is cut off at the Debye length pjj rather than extended to <», the ratio 
on the right side should estimate the effect if a typical v and the most 
Important level a' are used. 

This procedure overestimates the reduction of the width not only 

because it neglects the other perturbing levels a*, but also because the 

width vanishes even for a single level a' as soon as p approaches 

Pmin However, for v > v, p ^^^ is smaller than at v = v, 

so that a finite contribution to the width remains. If Debye shielding is 

lji 3 )ortant, that is, if the plasma frequency becomes of the order of the 

level spacing on between the initial (or final) state and the pertTorb- 
oa 

ing level a’, the effect of Debye shielding should be subtracted from the 
contribution of each perturbing level separately, and before the thermal 
average is taken: 


4 “ IX. »> 

A comparison has been made between w„ according to equations (lO) and (U) 
for the He I 58?6 line. Figure 3 shows the ratios w/N and s/w versus 
electron density. The estimated ratio w/H^ is also shown and can be seen 
to be reliable for width reductions of about 25 percent in this instance. 

The shift is considerably more sensitive to Debye shielding, to particular, 
Debye shielding must be considered if the first order independence of the 
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shift /width ratio is to be xitilized for electron temperatiire determinations, 
'even if this plasma thermometer, first suggested by D. D. Burgess (ref. 6), 
is calibrated experimentally. 

The measxrrements of axgon H broadening (refs. 3 and k) were carried 
out at 10^^ electrons/ cm^ density where Debye shielding is negligible for 
the He 5876 line which served for the density determination. 

IV. SUMMART 

In summary, the facts reported today are refinements to the mathematics 
involved in the GBKD theory and point to a need for line broadening calcu- 
lations which include Debye shielding. Such calculations are in progress at 
Langley. It is also concluded that the GBKO theory in its present form does 
fail for isolated lines of singly charged ions. The development of a more 
satisfactory theory is in progress. 
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